Damping law of photocount distribution in a dissipative channel 
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For a dissipative channel governed by the master equation of density operator dp/dt = 
k (2apa) — o) ap — pat a) , we find that photocount distribution formula at time t, p(n,t) = 

Tr{p(t) : ($a t a)"e" i ' ,, '/fi!: j , becomes p(n,t) = Tr \p(0) : (£ e - 2K Va)" e^ 2 ^ a " a /n\ : ] , as 

if the quantum efficiency £ of the detector becomes £e _2fct . This law greatly simplifies the theoretical 
study of photocount distribution for quantum optical field. 



In nature, syetems we concerned usually are sur- 
rounded by thermo reservoir. Decoherence is an impor- 
tant topic in the fields of quantum optics, quantum com- 
puting and quantum information processing. Decoher- 
ence of a quantum optical field can be judged by photo- 
counting its final state comparing with its initial state. 
In this paper we discuss time-dissipation of photocount 
distribution in a dissipative reservoir. We hope to know 
photon counting law for p(t), when the initial density op- 
erator p = p(t = 0) evolves with time in a photon-loss 
channel. The corresponding master equation is [l[ 



dp 
~db 



k (2apat 



a) ap - 



pa) a) 



(1) 



where k is the rate of dissipation, [a,a+] = 1. We shall 
deduce a new damping formula of photon counting after 
an optical quantum field passing through a photon-loss 
channel. 

Converting p 's master equations into pure state 's evo- 
lution equations. In order to solve equation Eq.©, in- 
stead of converting the operator equation into its c- 
number equation by the Wigner function approach or 
by P-representation approach, we introduce thermo en- 
tangled state representation |2| 



ra T -r*at + atat |00), (2) 



which is complete, to make up the pure state \p), 

oo 

\p)=p\I), |/}^|T = 0} = e atat |00) = ]T|nn), (3) 

71 = 

where a+ is a fictitious mode accompanying the real 
photon creation operator a), |5) is annihilated by a, 

[a, at] = 1, \nn) = {atat) n |00> /n\, [a, at] = 0. From 
Eq.© it is easy to see that 

a\I) = 5+1/}, 
ot|J) = a\I), 
(aU) n \I) = (tfa) n \I), (4) 



which indicates that the transformation relation between 
the real mode and the fictitious mode. Acting the both 
sides of Eq.© on the state \I) , we have 



dt 



\p) = n (2apat — at ap — pa) a) \I) 



(5) 



Using Eq. © and noting that p is defined in the real space 
which is commutative with the fictitious operators, we 
can convert Eq.© into 



—— \p) — k (2aa — a^a — a 1 a) \p) 



(6) 



Its formal solution is given by 

\p) = exp{nt (2aa — at a — at a)} \po) , (7) 

where \po) = pa\I). Thus by introducing a fictitious 
mode a+ to be a counterpart of the system mode, one 
can convert density operators' master equations into evo- 
lution equations of pure states, which leads to the formal 
solution ©. 

On the other hand, noticing that the commutative re- 
lation [(a^a + ata)/2 1 aa] = — aa, thus we can use the 
operator identity 



e \(A+*B) = e \A 



exp [c (l 



B/r] 



(8) 



which is valid for [A, B] — tB, to decompose the expo- 
nential operator in Eq.© as [3] 



-2nt( a ^ a + a ^ a —aa) —nt(a^a+a^a} Taa 



(9) 



where T = 1 — er 2Kt . On substituting Eq.© into Eq.©) 
yields 

\p (t)) = exp [— Kt (a+a + a+a)] exp [Taa] \po) 

= exp [-Kt {at a + at a)] ^ —a n a n p \I%10) 

71=0 U ' 

Noticing that the system mode is commutative with the 
fictitious operators, and the transition relation shown in 
Eq.©, we see that 
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a n po \I) = P0 a n \I) = Po o) n \I) 



(11) 



and 



cxp [-nta)a] a"poa tn 1-0 
= a n po(J n exp [— reta^a] | J) 
= a"p a t ™ exp [-/eto+o] | J) , (12) 

thus we have 

p(t)\I) = e- Ktaa Y—a n p a^ n e- Kt ° d \I} 



n=0 



_«t«Ta n t« e -«*» T <* |J) . (13) 

z — ' n' 



Comparing with the standard form of the operator-sum 
representation Q, 



p(t)= Y^MmPvMl 



(14) 



m— 



we see that the Kraus operator for this damping process 
is given by 




(15) 



and we can prove the unitrity X)m=o M) n M m = 1. 

Dam/ping formula of Photocount distribution in a dis- 
sipative channel. Now we discuss how does the photo- 
count distribution formula of optical fields evolve in the 
photon-loss channel? 

The quantum mechanical photon counting formula was 
first derived by Kelley and Kleiner As shown in Ref. 

the probability distribution p (n) of registering n 
photoelectrons in time interval At is given by 



P(n) = 1r{p: 



(fata) 



-£aJa , 



(16) 



where : : denotes normal ordering, £ oc Ar is called the 
quantum efficiency (a measure) of the detector, p (0) is 
the single-mode density operator of the light field con- 
cerned at initial time. According to the operator-sum 
representation at time t in Eq.(|14p. we have 



(17) 



p(n,t) = TrL(t):^^e-^: 



= Tr{p O(t)}, 
where we have set 



0(t)=^Ml.^Le-^-.M m . (18) 



Substituting Eq. (fT5)) into Eq. (IT8)) and using the following 
operator identities: : e^^ a : = exp [aJa In (1 — £)] = 

(l-£)° ta and e Aata at e - AQta = a t e A , e Aa1 a ae~ Aa1 a = 
ae~ A , we can put Eq. (fl"8|) into the following form 



6{t) 

^ T m + , t (&a) n e t , t 
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n! z — ' to! 

s_ e -2retn i a \n+m. e [ e [ 1 "< 1 -0- 2 " t ] _i] a t a . a "+"f]_g) 
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m=0 



Further using the property that a and are commute 
within the normal ordering symbol : :, we then remember 
T = 1 - e~ 2Kt to have 

6(t) 



= — -e a 
n! 



m=0 



ml 



-2kA 



a tn : exp{-^e _2,t *a t a}: a™ 



n! 



Substituting Ea.pD]) into (02) yields 



p (n, t) = Tr ^ p 



(^e- 2Kt ata)' 



(20) 



(21) 



Comparing Eq. (|2lj) with Eq. (jT6|) we obtain the law of 
damping formula of photocount distribution in a dissipa- 
tive channel, that is 



p(n,0) = Tr^(0) 



(£at a )" 



p(n,t) = Tr^p : 



^ ■ '— e -^ : Y22) 



m=0 



namely, this formula exhibits the law as if the quantum 
efficiency £ of the detector becomes £e~ 2K *. So from 
Eq. (f2~Tj) one can see that the photocount distribution is 
related to the decay of quantum efficiency £' = £e~ 2Kt . 
This is our main result. 

Example. Now we check the formula through an ex- 
ample. When the initial state is a pure number state, 
p (0) = | to) (to I , using the following operator identities 
a tm a ™ = N(N -1)---(N -m+1), where N = a^a, 
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then the photocount distribution is 



then using Eq. (fT6"| and Eg. ((23)) we have 



pM) 
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(m| a 1 n+l a n+l \m) 



^r- (m\ N (N - 1) • • • (N - n - I + 1) |m) 

n ' z=o 
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1=0 



(-0' r 



Z! n! (m — n — Z)! 



n! (to — n)\ 



(23) 



which indicates that the probability of a photon being 
counted during the period At is the quantum efficiency 
£, thus the probability of counting n out of to photons 
is proportional to the probability of counting n pho- 
tons £ n times the probability of counting m — n photons 
(1 — £) ™. Then according to the new formula Eq. ([22|) 
we immediately have the probability distribution of op- 
tical filed | to) after damping at time t, i.e., 



(f e - 2Kt ) n 



nl (m — n)! 



(24) 



Now we check it through another approach. We see that 
using (H~4l) - (p~5|) the initial state |m) (to| evolves into a 
binomial state, 



\m) (m\ — » e 



-IrnKt 



™ to! (e 2Kt - 1)' , 



1=0 



l\(m-l)\ 



(25) 



p{n,t) 

£^^<m-l|:^«^'- 



^ Z!(m-Z)!e 
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ml [e ZKZ - 1 



l\n\ (m — n — l)\ 



r (i - o 



m—l—n 



1=0 

n\ (to — n)\ ' \to — n — Z 



E( m ^,)(e-^i)'(i-e) 



rn —n—l 



n\ (to — n)\ 



(26) 



which is the same as Eq. (|24[) . Thus our new formula is 
correct. 

In summary, we conclude that for a dissipative system 
characteristic of amplitude damping, the photocount dis- 
tribution at any time can be considered as the one at 
initial time with the time evolution quantum efficiency 
£e~ 2tit , shown in Eq. (|2"2"|) . This law greatly simplifies the 
theoretical study of photocount distribution for quan- 
tum optical field. In addition, in Ref . [Toj . we have de- 
rived two new quantum-mechanical photocount formulas 
which build the relation from Wigner function and Q- 
function to photocount distribution. Thus our result can 
be applied directly to the case. 
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